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CN . We present the results of an analytical calculation of the distribution of gluons in the 
state immediately following a heavy ion collision in the quasi-classical limit of QCD given 
>■ I by McLerran-Venugopalan model. We show that the typical transverse energy Et of the 
^ ' produced gluons is of the order of the saturation scale of the nuclei Qs, as predicted by 
■ Mueller. We demonstrate that due to multiple rescatterings the obtained gluon distri- 
^ . bution remains finite (up to logarithms of k±) in the soft transverse momentum limit of 
k± <^ Qs- We compare the gluon spectrum in the single nuclear wave function before 
the collision to the spectrum of the produced gluons. The total number of produced glu- 
^ ' ons turns out to be proportional to the total number of gluons inside the nuclear wave 

function before the collision with the proportionality coefficient c ~ 2 In 2. 
Oh! 

r-| , The results presented here are based on reference |I| . Let us consider a collision of two 
ultrarelativistic nuclei characterized by the saturation scales Q^i and Q^2- We are going to 
write down an analytical expression for the classical distribution of produced gluons which 
resums all multiple rescatterings of the gluons with the nucleons in the colliding nuclei. 
Formally this corresponds to resumming powers of both Q'^i/kj_ and MM- The 

problem is similar to the problem of finding the classical gluon field of two colliding nuclei 
formulated in [Q, with the exception that here we are interested in the gluon production 
cross section. 

The classical field of two nuclei in the forward light cone has been found in the usual 
perturbation theory to order in |Q. The answer for the gluon multiplicity distribution 
is the first (lowest order) term of the expansion in powers of Q^i (see Eq. (^) below) 
and corresponds in this sense to proton-proton scattering. 

A gluon production cross section for a slightly more complicated case of proton-nucleus 
interactions was derived in 0. Formally that cross section includes all powers of Ql^/k\ 
keeping only the leading power of Ql2/^±- The calculations in ||^ were done in covariant 
gauge, where multiple final state rescatterings play crucial role in gluon production. 

In |]1| the same pA process was considered in the framework of light cone perturbation 
theory in light cone gauge. Knowing the answer for the gluon production cross section from 
allowed us to determine diagrams which reproduce this answer in the light cone gauge. 
We demonstrated that an entirely different set of initial state interactions is important 
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Figure 1. Diagrams contributing to the gluon production in nucleus- nucleus collisions in 
the = light cone gauge. 



there. Generalizing the pA results to the AA case we end up with the diagrams shown 
in Fig. m. Summation of the graphs of Fig. |l] yields us with the following multiphcity 
distribution of the produced gluons ^ 
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Eq. (|I|) is our main result. It provides us with the number of gluons produced in a heavy 
ion collision per unit transverse momentum phase space, per unit rapidity interval at the 
given impact parameter h yielding the initial conditions for possible thermalization of the 
gluonic system at later times. 

To explore the properties of the distribution (|I|) we first perform the transverse coor- 
dinate integration in the logarithmic approximation (see [jl|), which would transform it 
into the following simplified form 
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Figure 2. Distribution of the produced gluons given by Eq. (|^) for a central AA collision 
as a function of k/Qg with the transverse area of the nuclei being Ss_ = SOfml The 
approximation of Eq. @ is valid only for k± < Qs- 



The distribution of Eq. is plotted in Fig. |^ as a function of k/Qg for the case of 
two identical cylindrical nuclei with Qsi = Qs2 = Qs and with the cross sectional area 
S'_L = 7ri?2 ^ 50 fm^ ^ 1250 GeV"l Note that Eq. (|) is vahd only in the not very 
large transverse momentum region k±<Qs . 

As one can see, the distribution in Fig. ^ remains finite as k±/Qs — > 0. If one takes the 
k± <C Qs limit of Eq. then the distribution goes to a constant 

dN^^ k^ 

as — ^ 1. (3) 



d?kd?hdy as^ii^ Qg 

This means that the exact expression of Eq. (|^) may only have logarithmic divergences in 
the infrared limit. Therefore, finiteness of Eq. (|^) at small transverse momenta demon- 
strates that multiple rescatterings are the reason the hadronic and nuclear single particle 
inclusive production cross sections remain finite in the soft momentum region. 

Since the distribution of Eq. includes the lowest order in diagrams in it, one 
readily derives that in the k^/Qg oo limit the distribution falls off as @] 

dN^^ Ql.Ql^ k^ 

as — > oo. (4) 



d^kd^hdy agk Qs 

The average transverse energy per produced gluon given by Eq. (|^) for the case of two 
identical nuclei with Qsi = Qs2 = Qs is 

<Er>= Qs ^ l.OQQs. (5) 

That is, the typical transverse energy of the produced gluons is of the order of the sat- 
uration scale, as was previously conjectured by Mueller in The saturation scale for 
a large nucleus scales as Q^ ~ A^^^ with atomic number justifying the use of the small 
couphng expansion for large nuclei [@,|^. 
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Figure 3. The gluon distribution as a function of k^jQs in one of the nuclei before the 
coUision (thin hne) and the gluon distribution after the collision (thick line). 



Fig. ^ compares the transverse momentum distribution of gluons in the single nuclear 
wave function before the collision (thin line) with the distribution of produced gluons after 
the collision (thick line). One can see that the typical transverse momentum of gluons 
gets broadened due to multiple rescatterings in the collision. The total number of gluons 
produced in the collision can be found by integrating Eq. (^ over k^. The result yields 
for two identical nuclei 
dN^^ ^ CpQl In 2 
cPbdy ttsvr^ 

The total number of produced gluons turns out to be proportional to the total num- 
ber of gluons in the wave function of a single nucleus before the collision with the 
proportionality coefficient 

c = 2 In 2 ^ 1.39. (7) 

The obtained value for the "gluon liberation" coefficient c = 2 In 2 is close to one, as was 
originally suggested by Mueller 0. 
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